Universal decay cascade model for dynamic quantum dot initialization 
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Dynamic quantum dots can be formed by time-dependent electrostatic potentials in nanoelec- 
tronic devices, such as gate- or surface-acoustic-wave-driven electron pumps. Ability to control 
the number of captured electrons with high precision is required for applications in fundamental 
metrology and quantum information processing. In this work we propose and quantify a scheme 
to initialize quantum dots with a controllable number of electrons. It is based on the stochastic 
decrease in the electron number of a shrinking dynamic quantum dot and is described by a nuclear 
decay cascade model with "isotopes" being different charge states of the dot. Unlike the natural 
nuclei, the artificial confinement is time-dependent and tunable, so the probability distribution for 
the final "stable isotopes" depends on the external gate voltage. We derive an explicit fitting for- 
mula to extract the sequence of decay rate ratios from the measurements of averaged current in 
a periodically driven device. This provides a device-specific fingerprint which allows to compare 
different devices and architectures, and predict the upper limits of initialization accuracy from low 
precision measurements. 
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Single electron charging effects have attracted much in- 
terest since the proposal of single electronics [l| and the 
possibility to fabricate nanoscale structures. In partic- 
ular, quantum dots (QD) connected to leads have been 
a standard model system for many years to study sin- 
gle charges in so called artificial atoms. Dynamic QDs, 
which are repeatedly formed and manipulated by time- 
varying confining potentials appear in particular in struc- 
tures proposed to study quantum information 0^3 • One 
of the issues to be addressed is the decoupling of the 
QD from the environment and at the same time allow- 
ing the fast initialization with a controllable number of 
electrons. A suitable method has been demonstrated by 
Kataoka et al. J|] which uses pulses of surface acoustic 
waves (SAWs) to populate or depopulate a QD that is 
well isolated from electron reservoirs. While this mech- 
anism is not yet fully understood, a more conventional 
approach [§] is adiabatic decoupling of the QD from elec- 
tron reservoir, keeping the voltage on the QD fixed in a 
Coulomb blockade valley separating the discrete charge 
states. This strate^ is limited by (i) non-adiabatic ex- 
citation of the QD {9\ due to necessary finite decoupling 
rate, and (ii) experimental difficulties in tuning the lead- 
dot coupling to zero without disturbing the electrostatic 
potential on the QD, if. 

Here we propose and quantify an alternative scheme to 
achieve initialization by allowing non-equilibrium relax- 
ation (backtunneling) from a QD being raised energet- 
ically above the Fermi level during the decoupling pro- 
cess. This process is known to play a role in several 
types of non-adiabatic current generation devices 



14l |. We identify scale separation in integrated time- 



dependent electron escape rates between the subsequent 
charge states as a precondition for low-dispersion initial- 
ization. In our decay-cascade model a voltage parameter 



V shifts the hierarchy of the decay rates and thereby 
tunes the target number of electrons, rto, and the corre- 
sponding error ((n — no)^)- The minimal error is then 
given by the particular QD implementation and fixed, 
for instance, by QD geometry. We proceed to analyze 
the results of recent non-adiabatic current generation ex- 
periments ISMlTjl and find them to be promising realiza- 
tions of the proposed initialization scheme. Based on the 
decay cascade model predictions we propose an indirect 
way to measure the initialization accuracy in these de- 
vices by extracting the decay rate hierarchy fingerprint 
from low precision measurements of their current-voltage 
characteristics. Finally, possible strategies for accuracy 
improvement are discussed based on the ways charging 
energy, temperature and barrier design affect the electron 
escape rates. 

Decay cascade model. The initialization process is 
shown pictorially in Fig. [T^. The relaxation rate of elec- 
trons in the QD, Frc = {RG)^^ , is reduced at a char- 
acteristic speed (3 — |rRc/rRc|- (Here R and C are 
the resistance and the capacitance of the QD with re- 
spect to its environment, respectively.) Simultaneously, 
ip{t) grows more negative. Ejection of electrons at rate 
Fad ~ \^^lEc\ is required for the electron number distri- 
bution Pn[t) to stay close to instantaneous equilibrium 
{Ec = e^ /C is the charging energy and e is the electron 
charge). Such adiabatic following becomes impossible 
due to insufficient escape rate at times beyond the non- 
adiabatic crossover time to defined by FrcC^o) = Tad- 
Crucially for our scheme, if Fad can be made much larger 
than (3 then excitation from the Fermi sea [§] can be ig- 
nored while the QD remains loose enough at t > to to 
"forget" the adiabatic initial condition P„(to) and yield 
eventual accuracy dictated by the decay cascade mecha- 
nism. 
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FIG. 1. (color online) (a) Overview of the model, (b) Main 
results: the variance, (n — (n))'^ , and the average (inset), (n), 
of the number of captured electrons n, for Un = 1 , Si — 
and 5„ = S as functions of control voltage V (in units of 5). 
All results are calculated using N = 5 and Eq. Q , except for 
the thin dotted line that shows Eq. UHll with TV = 2. The 
thick dotted line traces the position and the value of variance 
as S is varied. 



troduce more specific physical assumptions that connect 
the model with potential experimental realizations. 

Conditions for accurate initialization. Consider an 
additional assumption (which can be relaxed later) that 
time-dependence of r„(i) is the same for all n, namely, 
r„(t)/r„_i(f;) = e'^" = const. In this case the final 
{t — >■ oo) distribution P„ depends only on the dimen- 
sionless integrals Xn = J^^ r„(i) dt and is given by 
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The solution ^ is rather unilluminating, but can be in- 
vestigated numerically (see below). For now let us focus 
on the limit of cascade timescale separation: 



> X„+i > Xn > Xn-l > 



(7) 



In this limit Eq. ([6]) simplifies to Qk-i,k = —^ and Qn^ — O 
for all n < fc — 1. This in turn means that C„ = 1 for all 
0<n<N. The solution becomes Pn{oo) = e"'^" and 



P„(oo) 



,-x„+i for < n < iV . 



(8) 



For t > to we write down a general kinetic equation 

dPn{t)/dt = ~Tn{t) Pn{t) + T n+l{t) F„+l(<) , (1) 

Pn{to)^5n,N. lim r„(t) = 0, (2) 

where r„ is the decay rate of the charge state with n 
electrons on the QD. The empty dot [n — 0) is an ab- 
sorbing state, Fq = 0, and the distribution is normal- 
ized, X]^o Pn = ^- N is the initial number of electrons 
and Fjv(io) is identified with Fad (see Fig. [1^). The use 
of Eqs. (IT])-([2]) to describe QD initialization assumes (i) 
randomization of the microstate corresponding to a given 
n on a time scale r ^ F,^^ so that Markov approximation 
is justified and the transition rates F„ are well-defined, 
(ii) no additional loading of electrons into the QD after 
to, and (iii) sharp initial condition, Eq. ^ (the latter 
is non-essential). The system of equations ([l} is rather 
general, and has also been used in the context of dynamic 
QD evolution [l^, [l^l (although not for the initialization 
stage). 

A general iterative solution to Eq. ([ij is 



Pnit) 



'n'^^i-)drY^^,{t')Pn+l{t')dt' , (3) 



to 



where P^+iit) = S{t — to)/TN+i{to) ensures fulfillment 
of the initial condition ©. Our strategy is, firstly, to 
identify the general properties of Eq. ^ that result in a 
well-defined final electron number, and, secondly, to in- 



Our model is justified for F7v(io) |rAr/Fjv|, which im- 
plies Xn ^ 1 (for smoothly decaying rates). Thus there 
will be such integer no < N that 



-'^rio + l > 1 > Xng and Xn„ + 1 ^ ^no- 



(9) 



Under Eq. ^ it is sufficient to consider only three prob- 
abilities to be non-zero: 
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We see that the probability distribution is dominated by 
Pno 1 if -'^no+i ^ 1 ^ii'i -^no ^ 1- The meaning of 
this condition is simple: the state with no -t- 1 electrons 
is unstable enough to have decayed while no is stable 

Time-scale separation expressed by Eq. ([7]) can be 
taken into account directly in the most general solution 
([3|), without requiring same time-dependence of F„(i)'s. 
The mathematics of this derivation can be summarized 
as follows: (a) observe that F„+iP„+i = -Jj I]„j>„-Pm 
exactly; (b) assume that all P,„(i)'s with m > n com- 
pared to Pn(t) itself vary on a time-scale much closer to 
to, so that Tn+iit')Pn+iit') oc S{to-t') in Eq. (c) 
solve the resulting P„(oo) = e"^" (1 - X]m=„+i ^m) to 
get 



P„(oo) 



-X, 



-Xrr 



fii) 



Finally, examining Ea. (|TT|) under condition (|9]) reveals 
difference from Eq. ([ro]) of at most X„„ which is reached 

at Xrin <C Xrin + l ~ 1. 
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Parametric control of decay rate hierarchy. Useful- 
ness of Eq. (HI) is limited unless dependence of X„'s on 
commonly accessible external control parameters can be 
established. Typically, escape rates depend exponentially 
on the height of the confining barrier, which in turn can 
be controlled by a gate voltage V. Tuning the decay rates 
^7i{t) by V would affect all decay rates simultaneously. 
Thus we propose 
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lnX„ = -a„y + V" 

^ ^l—l 



(12) 



Here a„ and (5„ are phenomcnological constants that can 
be readily extracted experimentally as discussed below. 

Figure [T}) shows the behavior of the first two moments 
for equal and fixed a„ ~ 1 and (5„ = (5 ((5i = is a mere 
shift of V) . For ^ > 6 the difference between the approx- 
imations (HI), (Uni) and (HU is negligible and Eq. is 
sufficient to describe the no-th step of a staircase (n) {V) 
regardless of N. The minimal variance {n — (rt))^ is 
achieved at optimal values of V that are easily found 
from Eqs. ^ and The minimal value of 1 - Pi (V) 
decays roughly exponentially with S. 

Possible experimental realizations. We expect to find 
good candidates for experimental realization of the 
proposed initialization mechanism among dynamic-QD- 
based electron p ump s since backtunneling relaxation has 
been found [lO|, [iJI to take place during certain parts 
of the pump cycle. The pumped current /pump consists 
of electrons captured from the source and subsequently 
ejected into the drain. It can be related to model predic- 
tion via 



/n 



ef (n) = e/ nP„(oo) 

^ — 



(13) 



(/ is the frequency which the pump cycle is repeated) 
provided that (a) ejection to the drain starts after the es- 
cape back to the source has stopped, and (b) the ejection 
is complete. In some experimental settings there is str ong 
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evidence that the latter condition can be ensured ^ ^ 

while in others 13, U | conditions (a) and (b) have been 
conjectured based on electrostatic modeling. 

We have used the ansatz (fT^ and the solution @ in 
Eq. (|13p to fit experimental data from various electron 
pump devices lB"17|, the results are shown in Fig. [2] In 
all cases q;„ were found to vary weakly with n, so we have 
set an = a constant for the entire voltage range for each 
device. Data set "a" in Fig. [2] was obtained in Ref. [15 
for a silicon nanowire metal-oxide-semiconductor field- 
effect transistor (FET) driven by voltage pulses, while the 
data labelled "b" correspond to a sinusoidal modulation 
of an AlGaAs/GaAs nanowire metal-semiconductor FET 
in the quantum Hall regime Yet another realization 
is exemplified by the data obtained in Ref. 17, for an 
AlGaAs/GaAs split-gate device ("c" in Fig. [2]), where the 
time dependent potential was generated by SAWs. The 
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FIG. 2. (color online) Upper panel, main graph: decay cas- 
cade model fits (solid line) compared to measured current 
7/ (e/) (pluses) from Fig. 2a of Ref. [isl ( a, magenta). Fig. la 
of Ref. llB fb. blue), and Fig. 1 of Ref.llTl (c, red), and to clas- 
sical simulation results from Fig. 6 of Ref. [ill (d, green). The 
graphs are shifted vertically for clarity. Insets are discussed 
in the text. Lower panels: fitted values of 5n for each case. 



potential and tunes the height of the confining barriers 
corresponding to the control voltage V . 

Figure [2] demonstrates particular advantage of a.n=a: 
(n — l)-th step in /pump(^^) has the length of 5„ on the 
scale of aV — 5\ . Thus the set of voltage- independent 
(5„ = XviXnl can be used a fingerprint of a par- 
ticular device, see the lower panels in Fig. [51 Plotting 
ln(— In (n — no -l- 1)) as a function of V serves as a quick 
test of the ansatz since the plateaux part dominated 
by X„(, must show up as a straight line. This is demon- 
strated in the lower right inset of Fig. [5] for the analytic 
(n) with b — ^ (same parameters as in Fig. [TJd) and in 
the upper left inset for the empirical data (same data sets 
as in the main panel). For contrast, we show by a green 
dash-dotted line an ad hoc fit with a sum of symmetric 
Fermi-like step functions; the deviation from the cascade 
model is notable. 

Our results for the second moment can be tested by 
measuring the low- frequency shot noise power spectral 
density [20] 



5 = 2e2/((n2)-(„)2) 



(14) 



Beyond the conditions (a) and (b) discussed above, 
Eq. ([H)) assumes [2^ that (c) the temporal width dr 
of the electron ejection current pulse is much less than 
the repetition period, dr/ <C 1. This regime has been 



voltage applied to the split gate superimposes the SAW- probed experimentally by Robinson and Talyanskii [21 



and we find good agreement with their data, see Fig. [3] 
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FIG. 3. Comparison to shot noise measurements on a SAW 
pump taken from Fig. 4a of Ref. [2ll . Upper panel: measured 
/pu,np/(e/) [♦ from Ipun,p{V), from Po{V) and P2{V)] and 
calculated (n) (solid line). Parameters fitted in Eq. Q with 
= 2 and ai = 02: ai = 92.04 V"\ Si = -240.7 and 
S2 — 6.252. Lower panel and central inset: experimental 
probabilities Pq (Q) and P2 (□) compared to Eq. Q (solid 
line) with parameters derived from the fit to Jpump- 



of the FET. Assuming that tunneling dominates in the 
latter case, (5„ is expected to scale proportional to the dif- 
ference in localization lengths between the ground states 
with n and n — 1 charges. 

In conclusion, we have proposed a scheme for accu- 
rate initialization of QDs and analyzed it quantitatively. 
We have shown that the model may be readily imple- 
mented using non-adiabatic pump architectures. It al- 
lows to extract a device-specific fingerprint which can be 
used to predict the results of a high precision measure- 
ment from a low precision characterization. In this way 
one can efficiently evaluate and optimize different pump 
architectures as required for fundamental metrology and 
adapt them for dynamic-QD-based quantum information 
devices. 
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Discussion and outlook. From the observations 
demonstrated in Figs. [2] and |3] we conclude that the de- 
cay cascade initialization scheme can be readily imple- 
mented experimentally. Several ways of controlling 
thus achieving more accurate initialization can be sug- 
gested: (a) tighter confinement to increase the charging 
energy, (b) employing more energy-selective barriers [22l | , 
or (c) lowering the local temperature. Suggestion (c) ap- 
plies when escape is determined by thermal activation. 
Note that this classical regime was investigated in elec- 
tron dynamics simulations of Robinson and Barnes [ll| , 
see trace "d" in Fig. [21 which demonstrates the universal- 
ity of the decay-cascade model. Within a classical inde- 
pendent electron picture, (5„ would be controlled by the 
difference in barrier height for the most energetic elec- 
tron, 6n ~ o,{En — En-i) whcrc En is the ionization en- 
ergy of a QD with n electrons, and a is an inverse effective 
temperature in this picture. The experiment of Ref. [15 
may have operated under these conditions since the re- 
ported temperature T — 20K is relatively high. The 
value of S2 = 6.4 from Fig. [2] then gives Ec = 11 meV in 
agreement the experimental estimate [l5| Ec = 10 meV. 
In other implementations lowering the temperature of the 
device may be not that effective: comparing SAW- [l3| 
and FET-based pumps jlj] (traces "c" and "b" in Fig. ^ 
respectively), we see a large difference in 62 of the first 
plateaux. Despite a similarity in QD area, temperature 
and material implementation the maximally achievable 
accuracy according to our model is {n — (rt))^ ~ 10^'^ 
for this particular SAW device versus 10^^ in the case 
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